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ABSTRACT 

We present some results from an ongoing investigation of the anisotropy of the 
two-point correlation function for optical galaxies. We have estimated £(7°j>,7T") 
from the Perseus-Pisces redshift survey, which is now virtually 100% complete for 
all morphological types to mg = 15.5. We detect strong distortions of the iso- 
correlation contours in £(r p , 7T"). These correspond to a global ('cluster' plus 'field') 

... _ 1 

pairwise velocity dispersion at 1 h Mpc, 0"(1), which is about a factor of two 
higher than the canonical value ~ 350 km s _1 obtained from the CfAl survey and 
recently from the 1.2 Jy IRAS redshift survey. However, the variance of this value 
is still rather strong even within the relatively large volumes surveyed (typically 
90° X 45° X 100 h -1 Mpc). From the largest volume explorable (126 h -1 Mpc 
depth), the best bet for the 'cosmic' value of 0"(1) seems to be between 600 and 700 
I km s _ . However, for a hopefully robust determination of its value we will have to 

wait for the analyses of the next generation of larger optical redshift surveys (such 
CO as, e.g., the now completed ESP 1 or the forthcoming Sloan Digital Sky Survey 2 ). 

o 

^ 1. Introduction 

Qh Peculiar velocities distort the maps of the galaxy distribution when radial veloc- 

O ities are used directly as a measure of their distance through the Hubble relation. 

£3 The observed distortions contain a wealth of important information on the statistical 

properties of large-scale motions, hence on the true mass distribution. In particular, 
the two-point correlation function in redshift space £(s) differs from that in real space 
£(r) in two respects: on small scales correlations are suppressed due to the virialized 
motions in rich clusters, which in redshift space elongate structures along the line 
of sight; on large scales coherent motions produced by infall into overdense regions 
or by outflow out of underdense regions have the opposite effect, i.e. they enhance 
correlations. The standard way to quantify these effects is to estimate correlations 
as a function of two variables, namely the two components of the separation vector 
on the plane of the sky and along the line of sight, indicated as £(r p ,Tv) . Projecting 
£(r p , 7r) onto the r p axis allows us to reconstruct the undistorted real space correlation 
function. In addition, a proper modeling of the observed £(r p ,Tv) allows us to charac- 
terize the pairwise velocity distribution function along the line of sight in terms of its 



moments. The first and second moments play a major role in describing the observed 
distortions on £(r p ,Tv) , being related, respectively, to the large-scale compression 
[the mean streaming ^12(7")] and to the small-scale stretching [the dispersion c(r)]. 
In general, both moments are a function of the separation r of the pairs. Here we 
shall report on new results concerning the second moment cr(r). 

The standard quoted value for cr(r) at small separations [in particular at 1 h -1 Mpc, 
cr(l) ] has long been that estimated by Davis & Peebles 3 on the CfAl survey, cr(l) = 
340 ± 40kms _1 . The recent determination of cr(l) from the IRAS 1.2 Jy redshift 
survey 4 seemed to provide a confirmation of such a value, with c(l) = 317^49 kms -1 . 
However, recent re-analyses of the CfAl survey 5,6 showed that the original value was 
rather sensitive to local corrections, like that for Virgocentric infall, pointing to a 
value about twice as high. As we shall show here, the analysis of an optically selected 
survey with a volume significantly larger than that of CfAl confirms that the value 
of cr(l) is most probably higher than the original estimate. The uncertainty of cr(r) 
due to variations in the sample size is however still quite strong, clearly indicating 
the need of larger surveys for obtaining a more robust estimate of its average 'cosmic' 
value. 
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Fig. 1. £(t°j>, 7r) for the sample with M < -19 and d < 79 h _1 . The heavy contour corresponds to 
^ = 1. For larger values of £, the contours are logarithmically spaced, with A log £ = 0.1 . Below, 
they are linearly spaced with A£ = 0.2 down to £ = 0. 

2. The Observed £{r p ,Tv) 

We defined r p and 7r as in 4 , and used the standard estimator 3 to obtain £(r p , tv) for 



volume-limited subsamples of the Perseus-Pisces (PP) redshift survey. (For a detailed 
description of the PP redshift catalogue, see e.g. 7 ). We stress the importance of using 
strictly volume-limited subsamples for surveys like PP (but also CfA2), where a 
dominant perpendicular structure is located near the peak of the selection function. 
In our case, this has the beneficial effect of suppressing the dominance of the Perseus- 
Pisces chain. It is interesting to note that in this case both the two-point correlations 8 
and the higher moments 7,9 agree between PP and CfA2. However, as we shall see, 
our results on velocity distortions indicate that these samples are still too small for 
providing the correct 'cosmic' mixture of cluster and non-cluster galaxies, necessary 
to estimate a global average value for cr(r) at small separations (similar to what can 
be obtained from large N-body simulations). 
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Fig. 2. Same as figure 1, but limiting the sample to RA < 3^, i.e. excluding the Perseus cluster. 

Figure 1 shows £(r p ,ir) for a sample limited to M < —19 (using H = 100 km s -1 
Mpc" 1 ) and d < 79 h" 1 Mpc, extending from 22 h to 4 h in RA and from 0° to 45° in 
DEC, containing 1159 galaxies. £(r p , it) has been estimated using a 2 h -1 Mpc binning, 
and for clarity the displayed maps have been smoothed with a 2-D Gaussian filter. 
Note the strong elongation of the contours along the line of sight at small r p 's, while 
on larger scales the opposite effect is evident. Part of this large-scale compression is 
most probably produced by streaming towards the PP supercluster. However, part 
of the effect is probably also due to the true spatial distribution of the pairs within 
the PP supercluster, not adequately compensated by other structures along different 
directions. In this situation, the two contributions cannot be disentangled using the 



redshift space data alone, thus discouraging any serious attempt to estimate the mean 
streaming on large scales. 

Figure 2 shows the effects on £(r p ,Tv) of simply restricting the sample to RA 
< 3 h 10 m , thereby excluding the richest cluster in the survey, the Perseus cluster. Al- 
though only 117 galaxies have been eliminated by this selection, the effect on £(r p ,Tv) 
is dramatic, significantly reducing the small-scale distortion. In the next section we 
shall quantify this change in terms of the implied pairwise velocity dispersion. 

3. The Small— Scale Pairwise Velocity Dispersion 

A proper modelling of the distortions induced on £(r p ,Tv) by peculiar velocities 
has to take into account the effects of, at least, the first two moments of the velocity 
distribution function /. Once a suitable functional dependence of / on the streaming 
^12 {f) and dispersion cr(r) is chosen, the problem reduces to finding a physically 
motivated parametrization of the scale dependence of these two quantities. This 
problem is discussed in detail in 4 . Here we present results obtained using the Davis 
& Peebles original streaming model, based on the similarity solution of the BBGKY 
equations. £(r p , 7r) is modelled as the convolution of the real space correlation function 
£(r) with the distribution function /. In the nonlinear regime, an appropriate choice 
for / is an exponential, so that one can write 



l + £(w) = #o [ + °° dy [1 + £(r)] - 7 J— exp|-V2tf 
J-°o V 2a(r) I 
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where r 2 = r 2 -\-y 2 , while the scale dependence of the streaming is given by the model 3 
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The other essential ingredient here is the assumption of a slowly varying cr(r), so 
that the second moment can be treated as a single free parameter. Given the observed 
£(r p ,Tv) , we first estimate the real-space correlation function by fitting a power-law 
model to the projected function w p {r p ) 4 . Then we fit the above model to slices 
through the observed £(r p ,Tv) at fixed r p , fitting simultaneously for the amplitude of 
the streaming F and the dispersion a . All the fits are performed using the maximum 
likelihood method discussed in 10 , based on 400 bootstrap resamplings. We limit the fit 
to small values of r p , where the velocity dispersion dominates over the mean streaming 
and the above model can be considered a reasonable approximation to reality. 

In Figure 3 we show the results of fitting the model to the first bin of the £(r p , tv) 
of Figure 1, i.e. the slice with < r p < 2 h -1 Mpc. The best fit values correspond to 
F = 0.54^^87 and to a dispersion a = 769^342 km s _1 . The fit is limited to ir < 10 
h -1 Mpc, i.e. where the power-law model for £(r) is a good representation of the 



data. However, the value of a is robust to changing the upper fitting limit to 20 or 
30 h -1 Mpc, with variations of a few tens of km s -1 . 

However, using the £(r p , 7r) map of Figure 2, i.e. excluding the Perseus cluster, the 
estimated a drops by nearly 200 km s -1 , as shown in Table 1. In the same table we 
list also the results obtained on two deeper sub-samples, up to 126 h -1 Mpc, again 
over the whole RA range. The aim here is to look for convergence of cr(l) as the 
volume increases. As can be seen, the estimated values still show strong fluctuations 
at the largest explorable sizes. 




Fig. 3. Fit of the model ^(t p ,7t) to the data from the M < -19 sample. Upper panel: Maximum 
Likelihood contours corresponding to the 68.3%, 95.4%, and 99.7% confidence levels on the two 
parameters taken separately. Lower panel: data points from the slice of the observed £(r p , 7T") , with 
the fitting model superimposed. 

Also, the lower 68% confidence limit is rather large for all the estimates but that 
from the M < — 19 "No Perseus" sample. Globally, one could dare to say that there 
is a tendency towards a value ~ 700 km s -1 , but this is certainly not statistically 
significant from these data. We are working towards improving the description of 
the streaming. However, as discussed by various authors 4 , the value of the second 
moment is rather robust with respect to variations in the streaming model, and we 
do not expect the best-fit values of cr(l) to change significantly. A full description of 
the work outlined here will be presented in a forthcoming paper. 

LG thanks J. Bahcall for hospitality at the IAS, where this work was initiated. 





Sample 


Ngal 


cr(l) (km sec *) 


M < 


-19.0 hd < 79 Mpc 


1159 


769t^ 




Same, no Perseus 


1042 


61311? 


M < 


-19.5 hd < 100 Mpc 


943 


927t^ 9 7 


M < 


-20.0 hd < 126 Mpc 


612 


71 7+190 
' -427 



Table 1. ML values for the pairwise velocity dispersion at separations between and 2 h Mpc 
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